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Tunable magnetic interactions in high-mobility 
nonmagnetic semiconductor heterostructures are 
centrally important to spin-based quantum tech- 
nologies. Conventionally, this requires incor- 
poration of "magnetic impurities" within the 
two-dimensional (2D) electron layer of the het- 
erostructures, which is achieved either by dop- 
ing with ferromagnetic atoms or by electro- 
staticallyprinting artificial atoms or quantum 
dots [H, El, 0, E]. Here we report experimen- 
tal evidence of a third, and intrinsic, source of 
localized spins in high- mobility GaAs/AlGaAs 
heterostructures, which are clearly observed in 
the limit of large setback distance (ss 80 nm) in 
modulation doping. Local nonequilibrium trans- 
port spectroscopy in these systems reveal exis- 
tence of multiple spins, which are located in a 
quasiregular manner in the 2D Fermi sea, and 
mutually interact at temperatures below 100 mil- 
liKelvin via the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) indirect exchange. The presence of such 
a spin-array, whose microscopic origin appears to 
be disorder-bound, simulates a 2D lattice-Kondo 
system with gate-tunable energy scales. 

Unintentional magnetic impurities are expected to 
be absent in high-quality nonmagnetic semiconductors, 
such as molecular beam epitaxy-grown GaAs/AlGaAs 
heterostructures. Contrary to this general belief, re- 
cent observation of a Kondo-like resonance in low- 
energy density-of-states of one dimensional (ID) quan- 
tum wires indicated existence of localized spin in meso- 
scopic GaAs/AlGaAs- based devices @|. Subsequently, 
evidence of localized spins was also reported in un- 
confined quasi-ballistic 2D systems 0, Hj], instigating 
the question whether localized spins arc intrinsic to 
GaAs/AlGaAs-bascd nonmagnetic semiconductors, and 
if so, what is the microscopic origin of these spins. A 
fascinating aspect of this problem is the possibility of a 
layer of mutually interacting spins: a system of consid- 
erable importance in studying various forms of magnetic 
ordering, quantum phase transitions and non-Fermi liq- 
uid effects (see Ref. [§] for a review). 

In thepast experiments with ballistic ID or 2D sys- 
tems [1,0,111, localized spins were detected with nonequi- 
librium transport spectroscopy. Pronounced struc- 
tures, commonly known as zero-bias anomaly (ZBA), 



in differential conductance (dl/dV) of mesoscopic de- 
vices close to zero source-to-drain electric potential 
(Vsd) were interpreted to be a consequence of spin- 
spin or spin-conduction electron exchange interaction. 
Here, we have augmented the nonequilibrium trans- 
port spectroscopy with perpendicular-field magnetore- 
sistance measurements, which not only confirm the ex- 
istence of multiple localized spins within high-mobility 
GaAs/AlGaAs heterostructures, but also reveals a strik- 
ing order in the spatial distribution of these spins that 
becomes visible over a narrow range of ri2D- 

Mesoscopic devices fabricated from Si monolayer- 
doped GaAs/AlGaAs heterostructures were used, where 
the 2D electron layer was formed 300 nm below the sur- 
face. A thick (w 80 nm) spacer layer of undoped AlGaAs 
between the dopants and the electrons provided a heavily 
compensated dopant layer with a filling factor / w 0.9. 
The resulting high electron mobility (~ 1— 3 x 10 6 cm 2 /V- 
s) provides a long as-grown elastic mean free path ^6 — 8 
/im, which acts as an upper limit to the device dimen- 
sions, ensuring quasi-ballistic transport. (Micrograph of 
a typical device is shown in Fig. la.) 

At low electron temperatures (T < 100 mK) and 
zero magnetic field, both equilibrium and nonequilibrium 
transport display rich structures well up to linear con- 
ductance G ~ 10 — 15 x (e 2 /h), as the voltage Vg on 
the gate is increased. For most devices the structures 
are strongest at carrier density U2D ~ 1 - 3 x 10 10 cto -2 
(Fig. lb), and consist of a repetitive sequence of two- 
types of resonances at the Fermi energy (Ep). This is 
illustrated in the surface plot of dl/dV in Figs, lc for 
device 1 of Fig. la. We denote the strong single-peak res- 
onance at Vsd = as ZBA-I, which splits intermittently 
to form a double-peaked ZBA with a gap at Ep, hence- 
forth referred to as ZBA-II. The illustrations of ZBA- 
I and ZBA-II in the inset of Fig. lc were recorded at 
points I and II in Fig lb, respectively. We define A 
as the half-width at half-depth of ZBA-II. While simi- 
lar nonequilibrium characteristics was observed in over 
50 mesoscopic devices from 5 different wafers, reducing 
the setback distance below ~ 60 — 80 nm was generally 
found to have detrimental effect on the clarity of the res- 
onance structures (for both ZBA-I and ZBA-II). often 
leading to broadening or complete suppression. 

The nature of suppression of such low-energy reso- 
nances with increasing T and at finite in-plane magnetic 
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field (~B||) indicate Kondo-like exchange (Figs, ld-lg) 
and hence presence of localized moments. A complete 
description can be obtained with the so-called "two im- 
purity" Kondo model, which embodies the interaction of 
an ensemble of localized spins within the sea of conduc- 
tion electrons 0, [II EE EI El- (See Ref. and 
Supplementary Information for arguments against alter- 
native explanation of the ZBA.) In the presence of anti- 
ferromagnetic coupling of individual spins to surround- 
ing conduction electrons, the zero-field splitting of the 
Kondo- resonance (ZBA-II) arises due to a Vc-dcpendcnt, 
oscillatory inter- impurity exchange J12, leading to the 
gap A ~ IJ12I at Ep. At certain intermediate values of 
Vg one obtains ZBA-I when |J 12 | < k B T. For ZBA- 
II, nonzero J12 results in a nonmonotonic suppression of 
dl/dV at low bias (|Vsd| ^5 A, as indeed observed exper- 
imentally (Figs, le and lg), while for ZBA-I this decrease 
is monotonic, and reflects suppression of single-impurity 
Kondo-rcsonance at individual noninteracting spins (see 
Figs. Id and If) [|. 

Mesoscopic devices showing clear resonances in the 
noncquilibrium transport also display a characteristic 
linear magnetoresistancc (MR) over the same range 
of U2U, when a small magnetic field (B±) is applied 
perpendicular to the plane of 2D electron layer. As 
shown for a different device, at high T (= 1.4 K), the MR 
consists of small peak-like structures at specific values of 
B±_ superposed on a parabolically increasing background 
(Fig. 2d), while at low T (= 30 mK), the magnetotrans- 
port breaks into quasi-periodic oscillations irrespective 
of the structure of the ZBA (Fig. 2a and 2b). We note 
that both these evidences indicate electron transport in 
a mesoscopic, quasi-regular array of antidots, which has 
been studied extensively in artificially fabricated antidot- 
lattices 16, 17, 18[. In the classical regime (high T), the 



peaks in MR correspond to commensurable cyclotron or- 
bits enclosing fixed number of antidots [la ], while at 
low T, quantum interference leads to the phase coherent 
oscillations in magnetoconductance, arising from trans- 
port along multiple connected Aharonov-Bohm rings as 
the inelastic scattering length exceeds sample dimensions 

The magnetotransport data shown in Fig. 2 allows an 
estimation of the inter-antidot distance R. In Fig. 2d, on 
subtracting the background, signature of commensurable 
orbits at B± w 0.08, 0.05 and 0.03 T corresponding to 
cyclotron radius of the electron encircling one, two and 
four antidots respectively (see inset). While this gives 
R ~ 500 nm (using the ri2p> ~ 1-31 x 10 10 cm~ 2 ), a more 
accurate estimate of R was obtained by Fourier trans- 
forming the low-T phase-coherent oscillations of Fig. 2a 
and 2b. In Fig. 2c, the power spectra calculated over the 
range to 0.065 T at both values of Vq show a strong 
peak at the frequency /ab ~ eR 2 /h as 105 ± 10 T _1 , 
corresponding to one flux quantum through unit cell of 
the antidot lattice (orbit b). This gives R ss 670 ±30 nm, 



which is consistent with the estimate from commensura- 
bility effect. Satellite peaks often appears in the power 
spectra, for example those at /ab — 190 T^ 1 (orbit c) 
and 50 T -1 (orbit a), which can be associated to specific 
stable orbits as indicated in the inset of Fig. 2a. R was 
found to be weakly device-dependent, varying between 
600 — 800 nm, but insensitive to the lithographic dimen- 
sions of the devices. 

Collectively, the observed nonequilibrium character- 
istics and low-field MR results indicate formation of a 
quasi-regular 2D spin-lattice embedded within the Fermi 
sea, where apart from the Kondo-coupling, the conduc- 
tion electrons would also undergo potential scattering at 
the lattice sites. The nature of such potential scatter- 
ing can be cotunneling [13], or scattering off the tun- 
nel barrier at the localized sites [l{|. In this framework, 
the RKKY exchange between the spins naturally leads 
to an oscillatory behavior of J12, where range function 
fy(2kpR) in the interaction magnitude reverses its sign 
with a periodicity of tt in 2kpR, kp = \f2im2D being the 
Fermi wave vector. Analytically [20(, 



A~ |J| ~ E F (Jp 2D ) 2 \V(2k F R)\ 



(1) 



where p2B is the 2D density of states and J is the ex- 
change coupling between an impurity spin and local con- 
duction electron. Fig. 3a shows the direct confirmation 
of this, where we have plotted A as a function of 2kpR 
for the device in Fig. lc. The clear periodicity of w n 
(within ±5%) in 2kpR can be immediately recognized as 
the so-called "2/cF-ff-oscillations" in the RKKY interac- 
tion, establishing the spin-lattice picture. 

The absolute magnitude of J12, and hence A, for a 
2D distribution of spins may differ widely from the sim- 
ple two-impurity RKKY interaction, and would be af- 
fected by frustrated magnetic ordering or spin glass freez- 
ing [2ll ] . as well as deviation from perfect periodicity in 
the spin arrangements [13]. Nevertheless, a framework 
for relative comparison of A in different samples can be 
obtained from Eq.[T]by normalizing A with Ep. As shown 
in Fig. 3b, adjusting for the experimental uncertainty 
in kp and R, A/Ep for four different devices with vari- 
ous lithographic dimensions can be made to collapse on 
the solid line proportional to modulus of pairwise RKKY 
range function over a wide range oi2kpR [20I ]. 

With known inter-spin distance, we shall now discuss 
two outstanding issues of this paper: (1) the microscopic 
origin of the uniform array of antidots, and subsequently, 
(2) emergence of the localized spins. A structurally in- 
trinsic origin of antidots in 2D Fermi sea of modulation- 
doped high-mobility GaAs/AlGaAs heterostructures can 
arise from long-range potential fluctuations in the con- 
duction band. Direct experimental imaging of disorder 
in similar systems reported typical distance between fluc- 
tuation to be ~ 0.5 — 1 /im, in excellent agreement to 
the magnitude of R in our devices [23[. In presence of 
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strong correlation in the dopant layer at large /, theo- 
retical investigations have also indicated a well-defined 
length scale in_the spatial distribution of the potential 
fluctuations 



24] 



A disorder-templated localized moment formation can 
then be envisaged through local depletion of electrons, 
analogous to moment formation in metal-semiconductor 
Schottky barriers We discuss this on the basis of 

three common representations of disorder and screening 
in high-mobi lity GaAs/AlGaAs systems, as schematized 
in Fig. 4a-c [25j] . At high n^o, the background disorder 
is linearly screened at all points, with local charge fluc- 
tuations |<5n 2 D| *C (ti2d) (Fig- 4a). With decreasing n2D, 
linear screening will break down locally at the maxima of 
slow potential fluctuations, where ri2P> becomes smaller 
than the local (rapid) density fluctuations of the dopants, 
resulting in the formation of single-particle localized 
states. Assuming a random distribution of the dopants, 
this is expected to occur at n% B ~ [(1 — ^ns/it] 1 / 2 
where rig = 2.5 x 10 12 cm~ 2 is the bare dopant density in 
our devices, and £ is the localization length. From /ab 
of orbit a in Fig. 2a we estimate £ ~ 150 nm, which gives 
?i^d ~ 1-9 x 10 10 cm -2 , which indeed marks the onset of 
strong ZBA (see Fig. lb). 

On further lowering of ri2Dj the system crosses over 
to the strongly localized regime (G -C e 2 //i), and the 
2D electron system disintegrates into puddles, which 
are often interconnected through quantum point con- 
tacts (Fig. 4c). While local spins can form at the 
point contacts 0, |2f|, the commensurability effect and 
phase-coherent magnetoconductance oscillations shown 
in Fig. 2 cannot be explained in such a picture, as they 
require extended and uninterrupted electron orbits. 

Our experiments thus outline a new microscopic 
mechanism of local moment formation in high-mobility 
GaAs/AlGaAs-based semiconductors with remote mod- 
ulation doping. This has serious implications on the pos- 
sibility of a gate-tunable static spontaneous spin polar- 
ization in mcsoscopic devices at low temperatures, and 
hinged on the Kondo-coupling of the localized moments 
to the surrounding conduction electrons. To verify such 
a coupling, we have estimated the ratio e/T for localized 
states using experimentally observed Kondo temperature 
T K , and that T K ~ (E F /k B ) exp(7re/2r) in the U -> oo 
limit, where e and T are the energy of single-electron state 
(with respect to Ep), and level broadening respectively 
(Fig. 4b), and U ~ e 2 /£ 3> Ep, is the on-site Coulomb 
repulsion. Taking the measured Ik ~ 265 mK at ZBA-I 
at point I in Fig. lc as an example, we find e/T ~ —2.1, 
which confirms the Kondo regime. 
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FIG. 1: Nonequilibrium characteristics: a, Picture of a set of two typical devices and electrical connections. Active area of 
a device is denned by the gate-covered region of the etched mesa. The data shown here was obtained from device 1, with 
device 2 and other side-gates kept grounded, b, Typical linear conductance, G, vs. gate voltage Vg (and electron density 
«2d) for device 1 at 30 mK and zero external magnetic field, c, Surface plot of differential conductance dl /dV of device 1 in 
Vsd — Vg plain. Each dl /dV — Vsu trace at a particular Vg was vertically shifted for leveling. The insets illustrate ZBA-I and 
ZBA-II-type resonances at points I and II in Fig. lb, respectively, d, Surface plot of ZBA-I in in-plane magnetic field {B\\). 
The dashed line shows a linear monotonic splitting with effective g-factor \g*\ « 0.5, which confirms the role of spin, e, Surface 
plot of ZBA-II in Bh. The single-impurity Kondo-behavior dominates above B\\ ~ A/g*fiB (the horizontal line), where A is 
the half-gap defined in text, f, Monotonic temperature (T) suppression of ZBA-I. g, T-dependence of dl jdV at ZBA-II. dl jdV 
is nonmonotonic in T for |Vsd| ^ A. 
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FIG. 2: Quantum and classical magnetotransport in perpendicular magnetic field (B±). Typical linear magnetoconduc- 
tance oscillation at a, a single-peak resonance (ZBA-I), and b, a double-peak resonance (ZBA-II). c, Power spectra of the 
magnetoconductance oscillations. The filled markers (blue: ZBA-I and red: ZBA-II) represent spectra obtained from the 
range \B±\ < 0.065 T, while the empty (blue) marker represents the spectrum (vertically scaled for clarity) from 0.065 T 
< B± < 0.15 T. The orbits corresponding to the peaks are indicated in the schematic, d, Four-probe linear magnetoresistance 
at 1.4 K for five electron densities from 1.22 (topmost trace) to 1.39 x 10 10 cm -2 (bottom trace), where well-defined ZBA's 
appear at low temperatures. Inset: Magnetoresistance after subtracting the parabolic background. The dashed lines, which 
denote various commensurate orbits, are computed using the average density of 1.31 x 10 10 cm -2 , and R m 500 nm. 



Siegert et al. 



6 




FIG. 3: RKKY indirect exchange and "2fcF_R-oscillations" : a, 2A from Fig. lc as a function 2k?R, where fcp is the Fermi wave 
vector, and R is the inter-spin distance obtained from magnetoconductance oscillations, b, 2A/Ep as a function of 2k?R for 
four different devices. Solid line is proportional to the pairwise RKKY range function (see text). Local disorder governs the 
experimentally attainable range of 2kpR within a given mesoscopic device. 
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FIG. 4: Schematic of background disorder and local moment formation. While the rapid fluctuations arise from dopant density 
fluctuations, the slow quasi-regular fluctuations indicates effect of strong correlation in the donor layer at large filling, a, large, 
b, intermediate, and c, low, carrier density regimes are shown. Our experimental results conform to the island-in-sea scenario 
of Fig. 4b. 



